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We explain simple laboratory experiments for making quantitative measurements of the Doppler 
effect from sources with acceleration. We analyze the spectra and clarify the conditions for the 
Doppler effect to be experimentally measurable, which turn out to be non-trivial when acceleration 
is involved. The experiments use sources with gravitational acceleration, in free fall and in motion 
as a pendulum, so that the results can be checked against fundamental physics principles. The 
experiments can be easily set up from "off the shelf" components only. The experiments are suitable 
for a wide range of students, including undergraduates not majoring in science or engineering. 
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I. INTRODUCTION 



The Doppler effect is a fundamental aspect of wave phenomena important in almost all areas of physics. The effect 
can be experienced in everyday life from fast moving cars or trains passing by and is also used as a practical method for 
measuring the velocity of objects such as baseballs and automobiles. It is also of crucial importance in astrophysics, 
where the receding or approaching velocity of an object is determined using the Doppler effect in light. The results 
are used as a basis for such fundamental aspects as Hubble's Law and the existence of black holes. We believe that 
the Doppler effect, given its importance and practicality, is a quite suitable theme as a basic physics experiment for 
all types of students. 

. — ■ In this work, we explain concretely experiments for quantitatively measuring the Doppler effect due to the motion 
of the wave source with gravitational acceleration. Our purpose is to design experiments for undergraduate students 
j^H *. including non-science majors. As such, we need an experiment that can be easily set up and safely conducted in an 
ordinary laboratory situation. Further, data acquisition needs to be robust and not sensitive to the environment. The 
!— students should be able to clearly understand what is going on and the results should be quantitative. Since we have 
a wide variety of students in mind, the analysis should not necessitate "higher" mathematics and desirably, the setup 
should use inexpensive components which are readily available. As explained below, we have been able to achieve 
all these goals. Since practical considerations are of import, we provide the details of the experiment, such as the 
| parameters and the equipment used. Due to its significance as a fundamental physics phenomenon, various forms of 
Doppler effect experiments have been proposed and conducted for some time: Doppler effect due to a wave source on 
a toy car—, on an air track 2 ^ 3 -^, in circular motion£&£, carried by a student^ and in free fal&i£ as well as the effect due 
I . to the motion of the receiver— have been studied. Demonstrations of Doppler effects from real cars have also been 
t—{ ' performe d 12 ' 13 . However, the majority of the past experiments seem either rather difficult to conduct, use ultrasonic 
ON . frequencies, not quantitative or require specialized equipment or circuitry. 

We explain the basic principle and the setup of the experiments in section ILT1 Section ITTT] contains the details and the 
| results from Doppler effect measurements of a freely falling source. We analyze the properties of Doppler effects from 
accelerating sources in section ITVl In the process, we clarify some non-trivial conditions and limitations in measuring 
Doppler effects. We explain the Doppler effect measurements of a pendulum source in section fVl 



II. EXPERIMENTAL PRINCIPLE AND SETUP 



In this work, we directly measure the Doppler shift due to the motion of the source, since this approach is most 
often used in various areas of fundamental physics and in practical measurements. The principle of the experiments is 
to measure the velocity of a moving wave source using the Doppler effect relation, f / f s = 1/(1 — v/c). The measured 
frequency is / s when the source is stationary and / when it is approaching the observer with velocity v. c is the speed 
of sound. From the time dependence of v, we can obtain the acceleration. While this principle can indeed be made to 
work, the simplicity is somewhat deceptive, even in principle. To measure a frequency at any given instant, a finite 
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amount of time T is necessary and the instantaneous velocity picture ceases to hold for large accelerations or for long 
measurement times, as explained in section HVl (T is different from the time of observation, during which we 

make multiple measurements of the frequency.) 

While the idea is simple, one main obstacle is that v is at most few% of c, due to practical considerations, including 
safety. Consequently, to measure v to 10% accuracy, for instance, we need 0.1% level accuracy in the measurement 
of /, which is not trivial. A theoretical limitation exists that to obtain a frequency resolution A/, time 1/A/ is 
necessary for the measurement. To obtain the desired precision, we use a relatively high frequency sound ~ 4 kHz 
and require A/ < 10 Hz. Consequently, the amount of time required for the measurement is T > 0.1s. A free fall 
from the height of 2 m takes 0.64 s, for instance, which is not that much larger than T. Increasing T to decrease A/ 
can conflict with measuring Doppler effects from accelerating sources, as explained in section ITVl T = Af /R so that 
A/ = 1/T = R/Af, where R is the sampling rate and Af is the number of samples. Any R can be used as long as R 
is higher than twice the frequency of interest. Af is chosen appropriately depending on the situation at hand. Since 
Af is fixed by T, larger / s can achieve higher relative accuracy. For the clarity of exposition to students and since 
generic components we utilize are not guaranteed to function properly outside the audible range, we use an audible 
sound source. All the experiments below were conducted with / s = 3.520 kHz, a clearly audible sound while not being 
too unpleasant. 

In the experiments, we need a simple way to measure frequencies in a short period of time with adequate accuracy. 
For this we use a microphone attached to a personal computer (PC) and a spectrum analyzer software. While most 
microphones work, one with high directivity facilitates the experiment. For spectral analysis, we used the GNU 
Octave software^ and an explicit example script for extracting spectra from a sound file is provided in Appendix 
An important technical point is to find an object that emits a sound at a distinct frequency steadily and a digital 
voice recorder is suitable for this task. It has the additional advantage that the sound can be synthesized on a PC, 
transferred to it, then played back. Synthesizing such a sound file can be accomplished by using GNU Octave^ or 
SoX 15 , for instance. A generic recorder with sufficient sound volume can perform the task and a means to directly 
transfer the sound from a PC is convenient. We also tested various simple buzzers for this purpose, but none of 
those tested could emit a sound with a steady enough frequency especially when in motion. Instead of a digital voice 
recorder, in some cases, a speaker attached to the same PC can be used to replay the synthesized sound, which can 
be somewhat cheaper. However, in this case, enough volume needs to be attainable and we need to use a software 
that does not emit the sound as it is recording, since it will otherwise interfere. Having an independent playback 
device such as a digital recorder, we find, is versatile, convenient and facilitates the experiment. The components 
we used are an inexpensive PC, Eee PC (ASUS, ~ $300), the microphone Audio Technica AT815b (~ $200) and 
the voice recorder, Olympus DS-51 (~ $100). All these components are generic and, in our experience, other similar 
components also suffice. 

III. A FREE FALL DOPPLER EXPERIMENT 
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FIG. 1: Examples of the measured sound spectrum, from the source at rest (solid) and when approaching the detector during 
free fall (dashed). The peak is Doppler shifted in the latter case. The subdominant peak corresponds to the reflection off the 
ceiling and is Doppler shifted as a receding source. The overall scale of |F(/)| 2 is arbitrary and Hann window is used. 

In this experiment, we let a sound source fall freely towards a microphone and measure the time dependence of the 
frequency. The experiment is similar in spirit to some previous experiments^^, though more quantitative. To obtain 
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the experimental results, we start recording on the PC and let the voice recorder drop (attached to a line to avoid 
it dropping all the way). The spectrum is analyzed at every At and the peak frequencies are obtained. Example 
spectra from a stationary and a freely falling source approaching the observer are shown in Fig. [T] Here, we used 
W = 2 15 , R = 192 kHz taking into consideration the constraints of section HVl Consequently, A/ = 6 Hz with an error 
Av = 0.6 m/s in v. An example of the time dependence of / is shown in Fig. [2] (left), where we can see the Doppler 
effect from a source falling repeatedly, after being pulled back up by the line. 
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FIG. 2: Doppler shifts measured from a freely falling sound source. (Left) Measured peak frequencies with At = 0.01 s. After 
4.2 s, the line is stretched to maximum and the source is pulled back up and falls again, which repeats. The frequency resolution 
A/ = 6 Hz is clearly visible. (Right) An example of an analysis performed in a student experiment. The same data as the 
previous plot, in the region t — 3.5 to t = 4.2 with At — 0.05 s. A linear fit to the Doppler shifted frequencies is shown and the 
slope corresponds to g = 9.4 ± 0.7 m/s 2 (ambient temperature 17° C). 



In the student experiments, we let the students read off the peak spectrum frequencies with At = 0.05 s or 0.1 s. Then 
the students plot these values against time and obtain the slope (Fig. [2] (right)). Using/ = f s /(l — at/c) ~ f s (l + at/c) 
(t = shifted here), they obtain the gravitational acceleration, g, from the slope. The students often discover the 
"bouncing effect" in Fig. [2] (left) during the analysis and it is educational for them to realize what they are observing. 
By using free fall for this experiment, the students can also clearly observe the constant acceleration due to gravity, 
which can contribute to the understanding of the gravitational force as well as the Doppler effect, g is measured 
within 10 % accuracy, enough to convince the student of the mechanism and the effectiveness of the Doppler shift 
measurements. The point here is confirm the understanding of the Doppler effect by obtaining g and its precision of 
g is not the main objective. Indeed, if measuring g were the primary objective, we can do so to 0.1 % accuracy using 
a simple pendulum. 



IV. DOPPLER EFFECT FROM ACCELERATING SOURCES 



In this section, we derive and explain the spectral properties of the Doppler effect from accelerating sources briefly. 
The practical problem is the following: Consider an accelerating monotonic source. We might expect the measured 
spectrum to have a peak Doppler shifted by the instantaneous velocity of the source, but a finite time is needed to 
make the measurement so no "instantaneous" velocity is measured. We provide quantitative criteria for observing a 
simple Doppler shifted spectrum and clarify what exactly can be measured. 

When a source approaching the observer with a constant acceleration, a, emits a signal at time t s , it is received 
by the observer at time t, and their relations are c(t — t s ) = L — at^/2. Here, without loss in generality, the initial 
distance between the source and the observer is L and the source is stationary when t s = 0. If the signal has a definite 
frequency / s , it can be reexpressed in terms of the observer time as 



sin(27r/ s t s ) = sin 



(1) 



Here, we assumed the source velocity v is such that v <C c and have kept the leading order terms. The signal is measured 
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by the observer from time t Q to t Q + T and its power spectrum can be obtained from the Fourier transformed signal, 

2 / , x 2 
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Here we defined a dimensionless parameter, 7 = a/(c/ s ). 7 is the relative change in the velocity in one period 
of the sound oscillation so that 7 < 1, We also defined G(±/) = C (x±,y ) , S(±f) = S(x±,y ) where x± = 
\f2rjf s {to — L/c) + (1 ± f/f B ) , y = \/2qf s T. These functions are related to the Fresnel integrals and their 

necessary properties are explained in Appendix [B] The power spectrum is 
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where we defined ip Q = 4(1 + / 2 // s 2 )/7- The spectrum is involved, but we first note that when / ~ / s , x+ — 
O 



-l/2\ 



^> 1. Therefore, using the asymptotics of the Fresnel functions in Eq. (|B3|) . G1/G0 = O (7 1 / 2 ) ,G2/Go = 
O (7) and both are much smaller than 1, when / ~ / s . We now analyze the dominant term Go and then show this 
suffices. As explained in Appendix [B] Go has a well defined peak when 



Vo = V^/fsT < Go «■ f s T 2 < ±C\ 



2a 



(3) 



where Go = 3.073490, numerically. y is the product of (the square root of) the acceleration relative to the sound 
velocity and T, both measured in the time scale l// s . Heuristically, this condition is reasonable; v is changing so that 
a simple peak does not exist when a is too large or T is too long. However, since both 7 and f s T are dimensionless, 
it is a priori not clear why this particular combination is the relevant criteria and this concrete computation enabled 
us to clarify the point. The condition puts an upper bound on T which competes with the precision A/ = 1/T, so 
that judicious balancing is needed. In all our experiments, T was chosen taking these considerations into account. 
The peak, when it exists, is at / pC ak//s — 1 = a/c(t — L/c + T/2). This Doppler shift corresponds to v for which 
the signal is received at the central time of the measurement, as expected. When y/2jf s T > Go, the spectrum has 
multiple peaks and the determination of the "peak" becomes ambiguous. 

We now analyze the effect of the subleading terms, Gi^, to / poa k- In general, when a function go(f) is peaked at 
/o, the peak for go(f) + e 5i(/) occurs at /o — eg'i(fo)/gn{fo)i to leading order in e. Therefore, we need to know the 
size of / dG\^{f) / df . Their individual terms are 0(7 -1 / 2 ) which can give rise to shifts comparable to the Doppler 
shift, so that a finer estimate of these terms is necessary. After some non-trivial computation, we obtain 
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The leading order terms in 7 cancel and f s dG\{f vca ^)/df — O (y Q ) , /srfG 2 (/pcak)/rff = O (7 1 ^ 2 ), using the asymptotic 
behaviors Eq. (|B3[) . Consequently, the contributions from G\ i2 to the relative shift in / pea k are 0{^y ), G(7 3 / 2 ), 
respectively. Since the relative shift in the frequency due to the Doppler effect is / pC ak//s — 1 > ^7/2 y /2, the 
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relative error in the shift due to Gi t % is at most C(7 1 / 2 ). We note that y a is at most of order one from Eq. ([3]) and 
typically not much smaller than 1 due to the required frequency resolution. 

We now illustrate these points in the free fall experiment in section Mil The condition for the existence of an 
unambiguous peak, Eq. ([3]), corresponds to f s T 2 < 166 s. In Fig. [3J we compare the cases when this condition is and 
is not satisfied both in simulations and in the analyses of the same experimental data. When condition Eq. (J3j is 
satisfied, a distinct peak exists in the spectrum. However, when the condition is not satisfied, a simple peak no longer 
exists even in theory and it is impossible to ascertain the peak frequency in the experimental situation. Condition 
Eq. (J3|) necessarily depends on the acceleration and for small acceleration or for constant speeds, we can use longer 
measurement times. For instance, if we use A/" = 2 17 = 131072 (which is not appropriate for measuring g), the error 
is A/ = 1Hz, corresponding to Av = O.lm/s. Errors Af,Av are inversely proportional to Af,T. However, due to 
Eq. (J3J, Af,Av can be decreased only as \/a. 
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FIG. 3: (Left) The spectra for R = 192 kHz, TV* = 2 J , f s T — 102.5 s, the simulated spectrum (top) and the experimental result 
(bottom). (Right) The spectra for R = 192 kHz, Af = 2 16 ,/ S T 2 = 410.1 s, the simulated spectrum (top) and the experimental 
result (bottom). The subdominant peak (left bottom) is due to reflection, similarly to Fig. [T] In all cases, no window function 
is used. 



Let us summarize the situation: When measuring constant velocity, we can reduce the error in v by using longer 
T. However, there will usually be practical limitations, such as one can not let the object move at constant velocity 
for an arbitrary amount of time nor can we pick up the sound for a long time. When a ^ 0, T needs to be small 
enough to satisfy Eq. ([3]) for the Doppler shift to be unambiguously observable. Further, if we want to measure o 
during time T \, s , we have an additional requirement that the Doppler shift needs to be larger than the frequency 
resolution. Letting ?7(< 1) be the relative precision with which we measure a, a7^bs/c > l/i]Coy / '2a/cf s , equivalent 
to Tq^s/T > 2(rjCoy ) • F° r larger a, 7ob s can be smaller, yet the attained velocity difference is larger. 



V. DOPPLER EFFECT MEASUREMENTS OF A PENDULUM SOURCE 



An experiment which provides more challenge in its analysis is the measurement of Doppler effects from a source 
used as a pendulum. We used the voice recorder strung by a line to a high beam as the pendulum and the configuration 
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FIG. 4: A top down view of the setup for a pendulum Doppler effect experiment. 



of the experiment is shown in Fig. |4j The strategy follows that used in the previous sections, exemplifying the ease 
with which we can adapt our methods to various Doppler shift measurements. The measurement is easier if relatively 
large velocities are attained. Consequently, we would like to a pendulum with a long length, £, and a large initial 
angle. Using typical values for a laboratory situation, £ — 2 m, #j n j t — vr/4, the maximum velocity is 3.5 m/s which is 
quite measurable, as seen in the previous sections. The experiment is not as simple as it may seem for the following 
reasons. First, if we use a large initial angle, the harmonic approximation to the oscillation is no longer quite valid. 
Furthermore, larger oscillations tend to die down rapidly. Smaller oscillations can be observed for a longer period of 
time, yet the velocities involved are smaller and the analysis becomes technically more delicate. 

Several approaches to analyzing this experiment are apparent. One simple approach is to take various time periods 
of the oscillation and compare them to theoretical predictions, as shown in Fig. [5] In these examples, £ = 2.32±0.01 m, 
the ambient temperature is 23°C. A lower sampling rate R = 48 kHz was used to reduce the data size and N — 2 14 
was used to satisfy Eq. The observed oscillation frequency agrees with the standard formula, \/g/l/ (27r). We can 
easily obtain g to 1 % accuracy or less and g — 9.8 ± 0.1 m/s 2 in this example. If the objective is solely to measure 
g, there is no distinct advantage over timing measurements of a simple pendulum. For our purposes, the pendulum 
source provides a variety of Doppler effects that can be quantitatively measured, including approaching, receding 
and accelerating sources all in one experiment, that can be checked against basic physics. Another advantage is that 
the time dependence of v can be obtained more globally than that in simple timing measurements. A precise global 
description of the data, which is not attempted here, needs to include the modeling of the decay in the oscillations and 
the anharmonic effects in the pendulum. As a Doppler effect measurement, we are quite satisfied that the oscillations 
can be seen clearly and agree so well with the theory. There are a number of aspects to this experiment and more 
intricate analysis can be performed if desired. For instance, the change in the period due to anharmonicity can be 
observed, as can be seen in the results in Fig. [5] This kind of experiment might be appropriate for more advanced 
students or for student projects. 



VI. DISCUSSION 



Considering the prevalence and the importance of the Doppler effect in basic sciences, along with its practical 
usefulness, we believe that easily realizable experiments in which the students can quantitatively measure the Doppler 
effect are highly desirable. In this work, we explained two such experiments, which are suitable for undergraduate 
students majoring in any field and require only generic, off the shelf components. In the experiments, the source 
undergoes gravitational acceleration so that the students can check the measured results against fundamental physics 
principles. The experiments are not difficult for students to perform and one can visibly confirm the induced shift in 
the peak of the audible sound. The setup we explained, as seen above, is quite flexible and the experiment can be 
conducted with the sound source attached to various objects in motion. It might be interesting to measure Doppler 
shifts of thrown objects, swung bats/rackets, and so on, using this method. 

In all our experiments, we picked out the highest peaks in the spectrum and no attempt was made to identify the 
"correct" peak, yet the Doppler effect could be measured quantitatively. The main reason for the non-ideal spectra and 
the occurrence of multiple peaks is the existence of reflections which inevitable occur in ordinary student laboratories. 
In our experience, objects that are approaching and therefore have higher peak frequencies tend to be easier objects 
for measurement, since the interference effects do not affect the highest frequency peak substantially. Our objective 
is to design student experiments which can be performed in ordinary laboratory situations and for such purposes, we 
designed the experiments with the data extraction that is robust and clear cut. 

The Doppler effect is a phenomenon that can be experienced in everyday life and is not hard for students to 
understand. One possible concern might be that students, particularly those not majoring in science or engineering, 
find spectral methods hard to comprehend. However, they are exposed to music and pitch in everyday life and most 
students understand the concept of frequency and harmonics, though not necessarily fully conscious of their spectral 
context. Therefore, the shift in the peak frequency caused by the Doppler effect is not a difficult concept for students 
to grasp. We believe, rather, that experiments with sounds allow students to pick up the idea of the spectrum and 
Fourier analysis intuitively and serves as an excellent introduction to spectral methods. In our experiments, we 
combine the spectral analysis of sounds from musical instruments with the free fall Doppler experiment in section HTT1 
in order for the students to first gain a good grasp of the relation between sound, harmonics and the frequency 
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FIG. 5: (Top left) Peak frequencies measured from a sound source swinging as a pendulum. (Top right) Peak frequencies and 
a function of the form (a + b e~ cot ) smuJo(t - to) with a = 11.5 Hz, b = 27.9 Hz, c = 0.045 s _1 ,o;o = 2.0125 Hz, to = —1.54s, 
which describes the data in this region well. The maximum velocity is 3.9 m/s. (Bottom left, right) Peak frequencies and a 
simple function, ai sincji(t — ii) with (oi, wi, t\) = (11 Hz, 2.04 Hz, 1.97 s), (7.5 Hz, 2.05 Hz, 2.34 s), respectively. The maximum 
velocity is 1.2, 0.8 m/s, respectively. A simple harmonic oscillation can be seen to describe the data in these time intervals well. 
fs = 3.52 kHz is also indicated in the plots. 

spectrum. Undergraduate students majoring in humanities and social sciences at Keio University currently perform 
these experiments in class smoothly and we are quite satisfied with the outcome. 



APPENDIX A: EXTRACTING SPECTRA FROM A SOUND FILE 



We provide an example script for GNU Octave (for Windows, Linux and other Unix variants) that extracts spectra 
from a sound file at predetermined time slices. This example script plots the spectrum at times 0, At, 2At, • • • , and 
displays the time (at the beginning of each measurement) and the peak frequency on the terminal, given a sound file. 
It can be modified to look at a particular region in time and window functions can also be incorporated. 



# ! /usr/bin/octave -qf 
[data,R] =wavread( ' f reef all . wav' ) ; 
dt=0.05; 
N=2~15; 

fmin=3400; fmax=3700; 
df=R/N; 

range=round( (fmin/df ) : (fmax/df)) ; 
f req=(range-l) *df ; 
printf("# t[s]\tfPeak[Hz]\n") ; 

for i=round ( (0 : floor (length (data) -N) /R/dt ) *R*dt ) ; 
spectrum=power(abs(fft(data(i+l:i+N))) ,2) ; 



# data, sampling rate from file 

# time between data points 

# number of samples 

# freq. plot range 

# frequency resolution 

# x axis 
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[specMax, jf]=max(spectrum) ; # peak frequency 

printf ("% . 5g\t% . 5g\n" , i/R, ( jf-1) *df ) ; # time and peak frequency 
semilogy(freq,spectrum(range)) ; # spectrum, semilog plot 
pause ; 
endf or 



APPENDIX B: SOME PROPERTIES RELATED TO FRESNEL INTEGRALS 



We define the following functions 

C(x,y)=C(x + y)-C(x) , S(x,y) = S(x + y) - S{x) . (Bl) 
Here C(x), S(x) are the Fresnel functions 16 , 

C(x) = J dt cos (|i 2 ) , S(x) = J dt sin (|t 2 ) . (B2) 

We need the asymptotic behavior of the Fresnel functions which are, to leading order.! 6 -, 



1 sin(7rx 2 /2) cos (nx 2 12] , ,. „, 1 cos(Vx 2 /2) sin (ttx 2 12] 

C(x ) = t; + — — — K -nr- L + ° ( x > s(x =- ^ — i-L L^_2 + o s-b). B 3 
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It is straightforward to find that the extrema of C{x,y) are at 



V 2n -y± y/8m-y 2 
2 + y ' 2 



and similarly for S(x,y), 



y , 2n -y±y/4(2m-l)-y 2 

x = , 

2 y ' 2 



(m > y 2 /8) (B4) 



(m>y 2 /8+l/2) . (B5) 



Here, m,n are integers. The function C(x,y) 2 + S(x,y) 2 has a maximum at x = — y/2 when y < Co = 3.073490 as 
seen in Fig. [6l where the value of Co was determined numerically. 
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FIG. 6: C(x,y) 2 (dashes), S(x,y) 2 (dots) and C(x,y) 2 + S(x,y) 2 (solid) for y = 1 (left top), y = 3.0735 (right top) and y = 4 
(left bottom). The maximum at a; = — y/2 disappears for y > Co. 



